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Abstract 

The hadronic shower longitudinal and lateral leakages and its 
effect on the pion response and energy resolution of ATLAS iron- 
scintillator barrel hadron prototype calorimeter with longitudinal 
tile configuration with a thickness of 9.4 nuclear interaction lengths 
have been investigated. The results are based on 100 GeV pion 
beam data at incidence angle = 10° at impact point Z in the 
range from —36 to 20 cm which were obtained during test beam 
period in May 1995 with setup equipped scintillator detector planes 
placed behind and back of the calorimeter. The fraction of the 
energy of 100 GeV pions at O = 10° leaking out at the back of 
this calorimeter amounts to 1.8% and agrees with the one for a 
conventional iron-scintillator calorimeter. Unexpected behaviour of 
the energy resolution as a function of leakage is observed: 6% lat- 
eral leakage lead to 18% improving of energy resolution in compare 
with the showers without leakage. The measured values of lon- 
gitudinal punchthrough probability (18 ± 1)% and (20 ± 1)% for 
two different hit definitions of leaking events agree with the earlier 
measurement for our calorimeter and with the one for a conven- 
tional iron-scintillator calorimeter with the same nuclear interaction 
length thickness respectively. Due to more soft cut for hit defini- 
tion in the leakage detectors the measured value of longitudinal 
punchthrough probability more corresponds to the calculated iron 
equivalent length Lp e = 158 cm. 



1 Introduction 



Due to limited dimensions of calorimeters one from important questions 
of calorimetry concerns the energy leakage and related with it the deterio- 
ration of energy resolution, appearance of tails in the energy distributions 
and ultimately the deterioration of the quality obtained physics informa- 
tion. In this article we report on the results of the experimental study 
of hadronic shower leakage effects on the pion response and energy reso- 
lution of ATLAS barrel hadron prototype calorimeter Because this 
calorimeter has innovative concept of longitudinal segmentation of active 
and passive layers (see Fig. |l|), the measurement of hadron showers leakage 
is of special interest 0. This investigation was performed on the basis of 
data from 100 GeV pion exposure of the prototype calorimeter at the H8 
beam of the CERN SPS at different Z impact points in the range from 
—36 to 20 cm with step 2 cm [Z scan) at incident angle = 10° which 
were obtained in May 1995. Earlier some results related with leakage for 
this calorimeter were obtained in ||, ||, 0. 

2 The Prototype Calorimeter 

The prototype calorimeter is composed of five sector modules, each span- 
ning 27r/64 in azimuth, 100 cm in the axial (Z) direction, 180 cm in the 
radial direction, and with a front face of 100 x 20 cm 2 0. The iron struc- 
ture of each module consists of 57 repeated "periods". Each period is 
18 mm thick and consists of four layers. The first and third layers are 
formed by large trapezoidal steel plates (master plates), 5 mm thick and 
spanning the full radial dimension of the module. In the second and fourth 
layers, smaller trapezoidal steel plates (spacer plates) and scintillator tiles 
alternate along the radial direction. The spacer plates are 4 mm thick and 
of 11 different sizes. Scintillator tiles are 3 mm thickness. The iron to 
scintillator ratio is 4.67:1 by volume. The calorimeter thickness at inci- 
dence angle 6 = 10° corresponds to 158 cm of iron equivalent (9.4 nuclear 
interaction length) ||. 

Radially oriented WLS fibres collect light from the tiles at both of their 
open edges and bring it to photo-multipliers (PMTs) at the periphery of 
the calorimeter. Each PMT views a specific group of tiles, through the cor- 
responding bundle of fibres. With this readout scheme three-dimensional 
segmentation is immediately obtained. 
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Tiles of 18 different shapes all have the same radial dimensions (10 cm). 
The prototype calorimeter is radially segmented into four depth segments 
by grouping fibres from different tiles. Proceeding outward in radius, the 
three smallest tiles, 1 -f - 3 are grouped into section 1, 4 7 into section 2, 
8-^12 into section 3 and 13 4- 18 into section 4. The readout cell width in 
Z direction is about 20 cm. 

Construction and performance of ATLAS iron-scintillator barrel had- 
ron prototype calorimeter is described elsewhere [[|, 0. 

3 Test Beam Layout 

The five modules have been positioned on a scanning table, able to allow 
high precision movements along any direction. Upstream of the calorime- 
ter, a trigger counter telescope was installed, defining a beam spot of 2 cm 
diameter. Two delay-line wire chambers, each with Z, Y readout, allowed 
to reconstruct the impact point of beam particles on the calorimeter face to 
better than ±1 mm 0. For the detection of the hadronic shower longitu- 
dinal and lateral leakages backward (80 x 80 cm 2 ) and side (40 x 115 cm 2 ) 
"muon walls" punchthrough detector were placed behind and at the right 
side of the calorimeter modules |§. Basic elements of u muon walls" are 
plastic scintillator detectors with dimensions 20 x40 x 2 cm 3 which are read- 
out by 2- inch photomultipliers EMI 9813KB. The tag of given (longitu- 
dinal or lateral) leakage is at least one hit in corresponding ll muon wall" . 
Due to the number of photoelectrons in any scintillator counter of walls is 
roughly 100 per minimum ionising particle ll muon walls" detected charged 
particles with high efficiency. As a result we have for each event 200 values 
of charges from PMT properly calibrated || with pedestal subtracted. 

4 Results 

30 runs contained 320 K events with various Z coordinates have been 
analysed. The treatment was carried out using program TILEMON 0. 

The scintillator detector planes behind and back of the calorimeter 
give us possibility to select the event samples at different conditions: "no 
leakage" , only "longitudinal leakage" , only "lateral leakage" , "longitudinal 
and lateral leakages" simultaneously. 

In this section the following issues are discussed: 
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1. punchthrough probability, 



2. energy leakage, 

3. the effect of leakage on energy resolution. 

First of all we determine the value of punchthrough probability. 

4.1 Longitudinal punchthrough probability 



By definition [§], |I0|| , [ |Tl|j the total hadronic punchthrough probability is 



the ratio of the number of events giving at least one hit in the punchthrough 
detector to the total number of trigger beam particles. It seems that the 
information needed is simple: hit or no hit. But there are some problems 
in definition of hit (see, for example, discussion JTD|). In Fig. ^| our ADC 
spectra one of u muon wall" counter (N- 8 in Fig. 3 ||) in u beam (top) 
and in 7r beam (bottom) are shown. Spectrum in ir beam look similar to 
simulated distribution for iron-scintillator calorimeter |Tl| as obtained by 
Monte Carlo calculations with GEANT (Fig. 14 from fH|). The region 
left from minimum ionising single particle distribution is related with the 
contribution of neutrons as punchthrough particles |10| . 
We used two cuts: 

1. ADCi > ADCj", where ADC^ — the beginning of Landau distribu- 
tion for z-counter, 

2. ADCi > (naturally after pedestal subtraction). 

Note that the results of cut 1 are not so much distinguished from a cut 
used in § ADC { > (< ADC >, t -3^). 

We think that cut 2 is more correct since it does not reject events with 
leakage. In following for the spectra analysis we use this cut. 

The results are given in Table [I], where in the last raw longitudinal 
punchthrough probability for different cuts corrected on value of accep- 
tance of the shower leakage detector (77 ± 4)% || are presented. 

As can be seen from this Table cut 2 is more soft relative to leakage and 
leads to decreasing of the events sample "no leakage" and to increasing of 
event sample with leakage. Especially the events sample with "longitudinal 
and lateral leakage" are increased (on 67%). 

Obtained value of punchthrough probability for cut 1 (18 ± 1)% agree 
with the one from ||. In the case of cut 2 obtained value (20 ±1)% more 
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Table 1: Percentage of the events and punchthrough probabilities for 
different types of leakages and cuts. 



Type 


Alias 


Cut 1 

% 


Cut 2 

% 


Cut 2/ Cut 1 - 1 

% 


no leak. 


nl 


72.0 


62.0 


-14. 


Ion. leak. 


11 


10.0 


9.4 


-6. 


lat. leak. 


lal 


14.0 


22.6 


61. 


Ion. & lat. leak. 


III 


3.6 


6.0 


67. 


all long. leak. 


lol 


13.6 


15.4 


13. 


punchthrough prob. 


PP 


18. ±1. 


20. ±1. 





correspond to calculated in f| iron equivalent length L Fe = 158 cm and 



the one for a conventional iron-scintillator calorimeter 12 



4.2 Energy response and leakage 

There are a few methods for evaluating of an energy leakage in calorimetry. 
For example, in []13| an additional "leakage" calorimeter was used for this 



purpose special. In [|14]], [[HJ the shower containment was measured by 
using the abundant longitudinal segmentation information. Since we do 
not have such possibilities the following method was used. We reconstruct 
the sum of initial energies of showers, E in , by using the detected energies 
of the event sample "no leakage", E nh and the fraction of these events, 
Nm/Nan: 

N all ]\J N „i 

E*4 = ^E^, (i) 

i=l n= i 

where N a u = N n i + Ni Q i, N n i — number of the event sample "no leakage", 
Ni i — number of the event sample "all longitudinal leakage" . 
The relative missing leakage energy is equal to: 



L 1 E% ggg j% N lol (< > - < Eg >) 

E£l ! El Ya=i E l in N all <E nl > 



(2) 



where En — energies of the event sample "longitudinal leakage" . 

In Fig.'s and f| two-dimensional spectra of energy responses as a 
function of Z coordinate and energy E are shown. Fig.'s |5| and [| show 
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the corresponding energy responses for events with all Z at different leak- 
age conditions. To map the energy in GeV scale the constant equal to 
100 GeV / <E nl > was used, where <E nl >= 514.2 pC is the mean energy 
response for event sample "no leakage". From these figures general be- 
haviour of energy response can be observed. It is seen that distributions 
for event samples "no leakage" and "lateral leakage" have almost Gaussian 
behaviour, the distribution for event sample "longitudinal leakage" have 
the clear low energy tail and in the distribution for event sample "lateral 
leakage" the maximum amplitude increases with increasing of Z. The ob- 
tained mean responses, relative resolutions as well as the values of leakages 
and tails are given in Table 0, where 

T _ < E n i > — < E,i > 

i = "no leakage", "longitudinal leakage", "lateral leakage", "longitudinal 
and lateral leakages", "all events". The estimate of tail is defined as an 
excess of the events over Gaussian curve in the region more than one 
sigma. 



Table 2: Responses, resolutions, leakages and tails for events with different 
Z in the range from —36 to 20 cm. 



Type 


% 


<E> 


u 

<E> 


<r n l 


L 


Low tail 


High tail 




Events 


GeV 


% 


% 


% 


% 


% 


no leak. 


62.0 


100. 


7.4 


0.0 


0.0 


0.0 


2.6±0.05 


Ion. leak. 


9.4 


91.0 


10.4 


41. 


9.0 


7.1±0.2 


0.0 


lat. leak. 


22.6 


96.8 


7.2 


-1.9 


3.2 


0.0 


1.1±0.05 


Ion. Sz lat. leak. 


6.0 


88.3 


10.7 


45. 


11.7 


6.1±0.2 


0.0 


all events 


100. 


97.7 


8.3 


13. 


2.3 


1.5±0.03 


1.3±0.02 



Fig.'s and [| show the energy distributions for event samples with 
various leakage conditions at Z = —8 cm. The characteristics of these 
distributions are given in Table |3|. The event samples — "any leakage", 
"longitudinal leakage", "longitudinal and lateral leakages" have the low 
energy tails, the event samples — "no leakage" and "lateral leakage" have 
the high energy tails. The events sample with leakage naturally have the 
low energy tail. The high energy tail in the event sample "no leakage" was 
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Table 3: Responses, resolutions, leakages and tails for the events with 
Z = —8 cm at various leakage conditions. 



Type 


% 


<E> 


(j 

<E> 




L 


Low tail 


High tail 




Events 


GeV 


% 


% 


% 


% 


% 


no leak. 


71.3 


100. 


7.3 


0.0 


0.0 


0.0 


2.7±0.2 


Ion. leak. 


11.1 


91.0 


9.9 


35. 


9.0 


6.7±0.7 


0.0 


lat. leak. 


14.2 


98.8 


7.1 


-4. 


1.2 


0.0 


1.7±0.3 


Ion. & lat. leak. 


3.3 


89.5 


8.8 


20. 


10.5 


12. ±2. 


0.0 


all events 


100 


98.4 


9.8 


7.9 


1.6 


1.5±0.1 


1.3±0.1 



explained in jOJ by contribution of showers with unusually large electro- 
magnetic component. The unexpected high energy tail in the event sample 
"lateral leakage" may be explained as these events are the events of type 
"no leakage" with some leakage unsufficient to cut the high energy tail. 

In Fig. ^| are shown the mean energy responses for events with differ- 
ent types of leakage obtained by averaging of energy spectra (top) and 
Gaussian fits (bottom) as a function of Z coordinate at different leakage 
conditions. In Table |] are given the results of averaging of these depen- 
dences in the uniformity ranges. 



Table 4: Responses and resolutions for the events at various leakage 
conditions. 



Type 


<E> 
GeV 


E G 
GeV 


RMS 

<E> 

% 


a 

E G 

% 


L 

% 


<Tnl 
% 


no leak. 


100. ±0.02 


99.7±0.02 


8.0±0.02 


7.4±0.02 


0.0 


0.0 


Ion. Ik. 


91.1±0.12 


94.0±0.08 


16.2±0.1 


9.7±0.07 


8.9±0.1 


31.0±0.8 


lat. Ik.* 


98.3±0.06 


98.1±0.06 


7.6±0.03 


7.3±0.03 


1.7±0.1 


-1.4 ±0.5 


all ev.* 


98.5±0.02 


98.8±0.02 


9.9±0.02 


8.0±0.02 


1.5±0.1 


8.1±0.4 



* For events with Z < 5 cm. 



The fraction of the energy leaking out from the backward side of this 
calorimeter calculated by the formula (0) amounts to (1.8 ± 0.03)% and 
agrees with the value 1.73% for L Fe = 158 cm measured in [|T5j|. 
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It should be noted that 15% of the events have the 9% energy longi- 
tudinal leakage and 1% of the events 50% of energy (« 50 GeV) leaking 
out at average. The latter estimate is extracted from the low energy tail 
in Fig. ^| (top). This fact must be taken into account in searching of new 
particles in future LHC experiments. 

We also considered the question concerning nonuniformity response of 
calorimeter. As can be seen in Fig. ^| the energy response as a function 
of Z coordinate from event sample "no leakage" is more uniform than the 
one for other event types. It is allows to estimate the more extended range 
of uniformity (from —36 cm to 20 cm) than in || which appears equal to 
0.9% (RMS). 



4.3 Influence of leakage on the energy resolution 

Fig. |TD| shows the relative energy resolutions obtained by Gaussian fitting 
of spectra (top) and the relative energy resolutions (RMS /<E>) obtained 
by averaging of spectra (bottom) as a function of Z coordinate at different 
leakage conditions. Fig. |ll] shows the same normalised to average value 
of (a /Eg) over the uniformity range for events without leakage. One can 
see that due to the tails the resolutions obtained by averaging are much 
greater (approximately in two times for events with longitudinal leakage) 
than ones obtained by Gaussian fitting. The results of averaged by Z in 
their uniformity range of Fig.'s |9], [K| |ll| are given in Table [|. As can be 
seen longitudinal energy leakage amounts 9%, but deterioration of energy 
resolution for the same case ct/Eg amounts 31%. The general degradation 
of the resolution with increasing of leakage is in agreement with earlier 
observations |13|], fl6f , fl7j , flSfl . Moreover, our energy resolution degrada- 



tion ( - cr " ° Vl '- > = 24% is in reasonable agreement with the parameterisation 



<y\x 



proposed by El on the basis of the data from CITF collaboration 



fo-*°) =0 . 9 /<fl>>-<fl> (4) 
a { , V < E > 

where < E >=< E nl > and a = o n i — energy and energy resolution 
for events without leakage, < Ei >=< En > and 07 = 07/ — energy and 
energy resolution for events with "longitudinal leakage". In our case for 
the value of energy resolution degradation from (|j) we obtain 27%. 

In the case of lateral leakage the unexpected inverse behaviour is ob- 
served: energy leakage leads to some improving of the resolution. Let us 
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consider this in more detail. In Fig. [12] two distributions of the lateral 
leakage are shown: lateral energy leakage (top), energy resolution (a /Eg) 
(bottom) for the event sample with lateral leakage as a function of Z co- 
ordinate. 

Fig. 13 presents the energy resolution as a function of lateral leakage 
for this event sample. As can be seen the energy resolution improves with 
increasing lateral energy leakage at least to the value of lateral energy 
leakage equal to 6% at Z = 18 cm where energy resolution is improving to 
18%. 

This phenomenon can be explained as follows. The hadronic shower 
consists of electromagnetic and pure hadronic parts and the electromag- 
netic part in lateral direction places in the central core fl7 l, So by 



cutting some lateral hadronic part we "improve" the shower properties, 
make it less fluctuating. However this may be the specific property of our 
calorimeter. 



5 Conclusions 

We have investigated the hadronic shower longitudinal and lateral leakages 
and its effect on the pion response and energy resolution on the basis of 
100 GeV pion beam data at incidence angle = 10° at impact points Z 
in the range from —36 to 20 cm. 

Some results are following: 

• The fraction of the energy of 100 GeV pions at G = 10° leaking out 
at the back of this calorimeter amounts to 1.8% and agrees with the 
one for a conventional iron-scintillator calorimeter. 

• Unexpected behaviour of the energy resolution as a function of leak- 
age is observed: 6% lateral leakage leads to 18% improving of energy 
resolution in compare to events with the showers without leakage. 

• The measured value of longitudinal punchthrough probability (20 ± 
1)% agrees with the one for a conventional iron-scintillator calorime- 
ter with the same nuclear interaction length thickness and with the 
earlier measurement ||. It also more correspond to calculated in || 
iron equivalent length Lp e = 158 cm. 
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Figure 2: Typical ADC spectrum of a ll muon walV counters in the \x 
beam (top) and in the ir beam (bottom) for counter N Q 8. 
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Figure 3: Two dimensional spectrum of energy response as a function of 
Z coordinate and energy E for various leakage conditions: all events (top), 
no leakage (bottom). 
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Figure 4: Two dimensional spectrum of energy response as a function 
of Z coordinate and energy E for various leakage conditions: longitudinal 
leakage (top), lateral leakage (bottom). 
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Figure 5: Energy responses for all Z at a different leakage conditions 
events (top), no leakage (bottom). 
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Figure 6: Energy responses for Z at a different leakage conditions: longi- 
tudinal leakage (top), lateral leakage (bottom). 



17 





Figure 7: Energy responses for all Z = —8 cm at a different leakage 
conditions: all events (top), no leakage (bottom). 
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Figure 8: Energy responses for all Z = —8 cm at a different leakage 
conditions: longitudinal leakage (top), lateral leakage (bottom). 
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Figure 9: Mean energy responses obtained by averaging of spectrum 
(top) and Gaussian fitting (bottom) as a function of Z coordinate at 
different leakage conditions: a) black square — no leakage, b) open square 
- longitudinal leakage, c) open triangle — lateral leakage, d) black circle 
— all events. 
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Figure 10: Energy resolutions (a /Eg) obtained by Gaussian fitting (top) 
and energy resolutions (RMS/ < E >) obtained by averaging of spectrum 
(bottom) as a function of Z coordinate at different leakage conditions: a) 
black square — no leakage, b) open square — longitudinal leakage, c) open 
triangle — lateral leakage, d) black circle — all events. 
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Figure 11: Normalised energy resolutions (a/E G ) obtained by Gaussian 
fitting (top) and normalised energy resolutions (RMS/ < E >) obtained 
by averaging of spectrum (bottom) as a function of Z coordinate at differ- 
ent leakage conditions: a) black square — no leakage, b) open square - 
longitudinal leakage, c) open triangle — lateral leakage, d) black circle - 
all events. 



22 




Z, mm 



Figure 12: The lateral leakage (top) and the energy resolution improving 
(bottom) for the events sample with lateral leakage as a function of Z 
coordinate. 
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Figure 13: The energy resolution improving for the events sample with 
lateral leakage as a function of lateral leakage for Z > — 5 cm. 
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